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ABSTRACT 

The first X-ray afterglow for a short (~ 30 ms), hard 7-ray burst was detected by Swift on 9 May 2005 
(GRB 050509b). No optical or radio counterpart was identified in follow-up observations. The tentative 
association of the GRB with a nearby giant elliptical galaxy at redshift z = 0.2248 would imply a total energy 
release E. 



7. ISO ' 



i 3 x 10 erg, and that the progenitor had traveled several tens of kpc from its point of origin, in 
agreement with expectations linking these events to the final merger of compact binaries driven by gravitational 
wave emission. We model the dynamical merger of such a system and the time-dependent evolution of the 
accretion tori thus created. The resulting energetics, variability, and expected durations are consistent with 
GRB 050509b originating from the tidal disruption of a neutron star by a stellar mass black hole, or of the 
merger of two neutron stars followed by prompt gravitational collapse of the massive remnant. We discuss how 
the available 7-ray and X-ray data provide a probe for the nature of the relativistic ejecta and the surrounding 
medium. 

Subject headings: binaries: close — stars: neutron — gamma-rays: bursts 



1. INTRODUCTION 

Classical 7-ray bursts (GRBs) naturally divide into two 
classes based on their duration and spectral properties 
dKouveiiotou et alJll993l) : short/hard (t < 2 s), and long/soft 
(f > 2 s) bursts. Through the impetus of the BeppoSAX satel- 
lite, it became clear that those of the long variety signal 
the catastrophic collapse of massive, rapidly rotating stars 
JWooslevl ri993) at high redshift (Metz ger et alJll997B . The 
nature of short events (about 1/3 of the total), is still unde- 
termined, but the merger of two com pact objects in a tight 
binary, as will occur i n PSR1913+16 Mulse & TavloJI 19751) 
and PSRJ0737-3039 JBurgav et alJ 120051) in 300 Mvr and 
85 Myr respectively, has long bee n considered a prime can - 
didate for a progenitor (Paczvhski 1986;: lEichler et all 1989). 
The short duration of the prompt 7-ray emission, however, 
precluded determination of accurate positions and follow-up 
observations, until now. 

A breakthrough came on 9 May 2005, when Swift suc- 
ceeded in promptly localizing GRB 050509b, a short burst 
lastin g only r 5 o ~ 30 ms JGehrels et alJl2005t IB loom etal] 
120051) . The fast response allowed for an accurate position de- 
termination and a rapidly fading X-ray source was located, 
falling below detection within ~ 300 s. For the next few 
days several multiwavelength observations were made, but 
unfortunately no optical or radio afterglow was detected. Al- 
though the issue of a host and its implications for the distance 
scale remain to be resolved, initial reports of a giant ellipti- 
cal galaxy at redshift z = 0.2248, lying only 10 away from 
the burst position 1 appear consistent with model expectations 
of compact binary mergers. This is because a compact ob- 
ject binary could take hundreds of millions of years to spiral 
together, and could by then — if given a substantial kick ve- 
locity upon formation — have traveled sev eral tens of k ilopar- 
secs away from its p oint of origin (see e.g., Bloom et a l J 19991 
llvanova et alJl2003l for population synthesis estimates). The 
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detection of GRB 050509b thus presents us with the unique 
opportunity, to which this Letter is devoted, to constrain this 
scenario, both from the prompt 7-ray emission and the after- 
glow. In §|2] we address the energetics and timescales which 
can be expected for the merger of two compact objects based 
on recent calculations, and compare them with the data for 
GRB 050509b. In §[3]we constrain the properties of the ejecta 
and the external medium by using the information available 
to us from both the afterglow and prompt emission, consider- 
ing both the distance scale of the tentative host galaxy and a 
higher redshift. Our findings are summarized in §@] 

2. ENERGETICS AND INTRINSIC TIME SCALES OF THE TRIGGER 



It has long been assumed (Lattimer & Schramm 1974J that 
the merger of a black hole-neutron star (BH-NS) or double 
neutron star (NS-NS) binary would result in the formation 
of an accretion disk with enough mass and internal energy 
to account for the energetics of a typical GRB, through tidal 
disruption of the neutron star in the former, or post-merger 
collapse of the central core in the latter. Calculations sup- 
porting this view have been carried out in the Newtonian 
regime, resulting in disks with « 0.3M Q , k T sa 10 MeV, 
p ss 1 n g cm" 3 , which could power a GRB (| R uffert et alJ 
Il996t iKluzniak & Lelfl998l iRosswog et alJl2002l) . General 
Relativity (GR) is certain to play a role, but gauging its ef- 
fects is not an easy task. The star could plunge directly 
into the black hole and be accreted whole in a matter of 
a millisecond (Miller 2005), precluding the production of a 
GRB 10 to 100 time s longer. Pseudo-Newtonian simula- 
tions ( Rosswog 2005) and post-Newtonian orbital evolution 
estimates (Prakash et al. 2004), however, reveal that the star 
is frequently distorted enough by tidal forces that disk-like 
structures and long, partially unbound tidal tails can form. 
The outcome is sensitive to the mass ratio q = Mns/Mbh 
and it appears that ro tating BHs favor the creation of disks 
jTaniguchi et alJl2005l) . For mass ratios q ~ 0.25 it is possi- 
ble to form a disk 2 , although of lower mass than previously 

2 For the 18 galactic BH binaries, an absolute lower bound is Mbh > 
'i.IMp), and for 8 of them (44%) , average values yield 6.5 < M-qh/Mq < 7.5 
IMcClintock & Remi llard 2004). 
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TABLE 1 

Disk formation in BH-NS mergers. 
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0.3 
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0.03 
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0.05 
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0.01 
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0.04 
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2.0 


0.03 
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0.1 


2.0 




0.02 



thought, ra d « 10" 2 M Q . 

To better estimate the mass of the disk (which will crucially 
affect the energetics) and the circumstances under which it 
may form, we have extended our study of merging BH-NS 
pairs using a pseudo-Newto nian potential in three dimen- 
sions (Lee & Kluzniak 1999) and summarize our new results 
in Table ^ A relatively narrow, but not unlikely range of 
parameters allows for the formation of a small disk, with 
rad w 3 x 10~ 2 M Q . Mass ratios higher that 1/3 are unlikely 
to occur, and if q < 0.1 only a wide, relatively cold arc-like 
structure is formed. The densities and temperatures in the re- 
sulting disks are p w 10 10 - 10 n g cm" 3 and kT w 2-5 MeV. 
We have considered the stiffness of the nuclear equation of 
state as a parameter by using polytropes with various indices 
in the range 5/3 < T < 2. The standard mass for the neutron 
star is 1 .4M Q . 

In the case of merging neutron stars (IShibata et alJ l2005). 
a low-mass disk (with « 1% of the total mass) may survive 
once the supra-massive remnant collapses because of gravi- 
tational wave emission on a time scale shorter than w 100 ms, 
and release up to 10 50 erg in neutrinos. In addition, the merger 
process and the collapse itself would likely p roduce a signal 
of their own (Rosswoa & Ramirez-Ruiz 2002). 

Once a disk is formed, the energy output depends on its ini- 
tial mass, mn an d temperature. We have recently calculated 
l lLee et alJl2005l) a realistic set of time-dependent models for 
their dynamical evolution, covering the typical duration time 
scales of short GRBs (up to one second). These 2D models 
make use of a Smooth Particle Hydrodynamics (SPH) code 
in azimuthal symmetry and include an accurate equation of 
state which considers photodisintegration, neutronization, and 
a relativistic Fermi gas of arbitrary degeneracy, as well as neu- 
trino cooling and finite optical depths to neutrinos. From the 
resulting neutrino luminosities we have computed the total en- 
ergy deposition that could drive a relativistic outflow through 
vV annihilation, assu ming a 1% efficiency at L v = 10 53 erg s" 1 
( Pophametal . 1999) and its duration. The results for various 
disk masses and effective a-disk viscosities are shown in Fig- 
ure [D(j omec l square symbols), along with the energy-duration 
curve for GRB 050509b constrained by the redshift. The to- 
tal output Euu ~ 10 49 [m d /0.03M o ] 2 erg, is roughly indepen- 
dent of the inferred duration, which increases with decreas- 
ing disk viscos ity since the overall evolution is slower (see 
iLee et all2005l) . The strong dependence on disk mass reflects 
the sensitivity of the neutrino emission rates on temperature 
(e^ capture on free nucleons dominates the cooling rate, with 
emissivity q oc pT 6 ). 

Magnetically dominated outflows may alternatively tap the 
disk energy through the Blandford-Znajek mechanism. Our 
estimates are shown in Figure ^ (round symbols), assuming 
equipartition of the magnetic field energy density and the in- 




time / ms 

FI G. 1. — Top: H istogram of observed short-hard burst durations taken 
from Paciesa s et al. 1 119991) . Bottom: Comparison of the energy-duration 
relation as a function of redshift for GRB 050509b (black line) with estimates 
from compact binary mergers. The square (round) joined symbols show the 
total isotropic energy release (assuming collimation of the outflow into = 
47r/10) and duration (fso) for up annihilation (Blandford-Znajek)-powered 
bursts, as computed from our 2D disk evolution models. The range in initial 
disk mass covers one order of magnitude and the effective disk viscosities 
are a = 10 _1 , 10~ 2 , 10~ 3 (left to right). Many of the estimates are lower limits 
because at the end of our calculations not enough mass had drained from 
the disk for the luminosity to drop appreciably. The stars correspon d to itv— 
driven outflows in NS-NS mergers, computed by Rosswoa & Liebendorfei 



ternal energy of the fluid in the inner disk. The energy flux 
is sensitive primarily to the equatorial flow density. It is thus 
initially roughly constant, then drops on an accretion (i.e. vis- 
cous) timescale. This explains the energy-duration correla- 
tion in Figure ^ Since the observed flux sets the thresh- 
old for burst detection, neutrino powered events will enhance 
the relative importance of shorter events (since E v p ~ const, 
then L„ p oc f -1 ), while magnetically dominated short GRBs 
more truthfully reflect the underlying intrinsic distribution 
(since Lbz ~ const, then E v ^ oc t). Relaxing the assump- 
tion of full equipartition will lower the total energy budget 
accordingly. The dependence on disk mass is different, with 
E BZ ~ 5 x 10 50 [TO d /0.03M Q ][a/10- 1 ]-°- 55 erg. Our estimates 
assume that whatever seed field was present has been am- 
plified to the correspondingly high values extremely rapidly. 
Whether this will actually occur is unclear, particularly for the 
shortest events, as the field can grow only in a time scale asso- 
ciated with proto-neutron star-like convection or differential 
rotation in the case of the MRI. 

3. CONSTRAINTS ON THE PROPERTIES OF THE EJECTA AND THE 
EXTERNAL MEDIUM 

The afterglow of GRB 050509b was detected by the Swift 
XRT during an observation which started 62 s after the burst 
and lasted 1.6 ks (Bloom et al. 2005) with a flux of Fx ~ 
7 x 10" 13 erg cm" 2 s" 1 in the 0.2-10 keV range at t = 200 s, 
and a temporal decay index awl .3^ 3, where F v oc t~ a . The 
numerous upper limits in the optical and few upper limits in 
the radio are not very constraining for the theoretical models 
(see Bloom et al. 2005). The fact that the X-ray flux was 
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already decaying at t > 60 s implies fd ec < 60 s, where 



fdec = (l+z)^=42(l+z) 
ZCL 



Esi 
n 



1/3 



100 



-8/3 



(1) 



and /?dec = (3-£ , k,iso/47rnm / ,c 2 ry) 1 ' /3 are the observed time and 
radius where the outflow decelerates significantly, To is the 
initial Lorentz factor, n = «o cm" 3 is the external density and 
£ic,iso = 10 51 £5i erg is the isotropic equivalent kinetic energy. 
That is To = 87[f dec /(l +z)60s]- 3 / 8 (£ 5 i/«o) 1/8 . 

3.1. Prompt Emission from Internal Shocks 

Internal shocks typically occur at a radius /?is « 2Tlct v 
where t v is the variability time. Since GRB 050509b had a sin- 
gle peaked light curve (Gehrels et al.l2005l) . t v ~ Tgrb/CI +z) 
where Tqrb ~ 30 ms is the observed burst duration. The 
Thompson optical depth is tj = E^ aj /4TTR 2 m p c 2 To. To see 
the prompt emission we need tt(Ris) < 1 implying 



r > ioo^ 5 



30 ms 



-2/5 



(2) 



For internal shocks, the vF v spectrum peaks at 



: hV m 



1.3g 2 



t B,-2 t e, 



r 



f7l/2 (30ms) 3 / 2 / 

451 tT W ^100 



keV 



(3) 

where eg = 10 2 €b,-2 and e e = 0.1e c .-i are the fractions of the 
internal energy behind the shock in the magnetic field and in 
relativistic electrons, respectively, g = 3(p—2)/(p- 1) and p is 
the power-law index for the electron energy distribution. Eqs. 
|2]and[3]imply 




e 1/2 e 2 E 

^.-2^-1^51 



1/10 



(—) 
V30ms/ 



30 ms/ \30ms 



-1/2 



keV 



(4) 



The S wift BAT spectrum is vF „ cx i/°~ 5±0A in the 15-350keV 
range ( Barthel mv et al 120 05). implying E p > 300 keV, which 
is hard to achieve for internal shocks (see Eq. |4}. Possible 
ways of increasing E p are if (i) the internal shocks are highly 
relativistic, rather than mildly relativistic as assumed above, 
or (ii) only a small fractio n of the electrons are accelerated 
to relativistic energies (e.g. Ramirez-Ruiz & Llovd-Ronnine 
120021) . It is not clear how likely either of these options is. The 
constraints on the physical parameters in the internal shocks 
model are summarized in Fig. [2] 

3.2. Prompt Emission from the External Shock 
In this case fdec ~ Tgrb ~ 30 ms, implying a very high To « 
1500(£5i/no)^ 8 - Also, E p = max(hi> m ,hi> c ) where 



hv„, 



i^e^e^^s^fdec/SOms)- 3 / 2 MeV , 



(5) 



/! i / e = 25(l+Fr 2 e B 3 / 2 no 1 £ 5 i /2 (fdec/30msr 1/z keV , (6) 

and Y is the Compton y-parameter. The value of E p is reason- 
able and independent of n for v c < v m . This requires, however, 
sufficiently high values of n and Esso- 
in the external shock model the prompt emission and the 
afterglow are produced in the same physical region. It is thus 
instructive to check whether the extrapolation of the flux in 
the prompt emission to the XRT observation at t i=s 200 s re- 
produces the observed flux. The prompt fluence was / « 
(2.3 ± 0.9) x IP" 8 erg cm~ 2 in the 15-350keV BAT range 
iBarthelmv et al.l2005l) implying a 7-ray flux of f 7 (20ms) « 



z = 0.22 




/n=ioo f 



n ext /crTr 

FIG. 2. — Constraints on the isotropic equivalent kinetic energy, E^ ls0 , 
and the external density, n, for the redshift of the tentative host galaxy (z = 
0.2248; upper panel) and for z = 3 (lower panel). Dashed lines labeled by the 
value of the initial Lorentz factor, Vq, bound the regions of allowed parameter 
space (in the direction of the arrows). These limits apply only if the prompt 
emission is from internal shocks, and are derived from the requirements that 
'dec < 60 s and tt(Ris) < 1- The shaded region is that allowed from the X- 
ray flux at the Swift XRT observation (f rj 200 s) for a reasonable range of 
values for the micro-physical parameters: 2.2 < p < 2.5, 0.03 < ee < 0.03, 
10 -3 < eg < 0.1 (this is independent of the model for the prompt emission). 
The 'plus' symbols show the location of the exemplary models given in Table 
3 of Bloom et al. (2005). 



10 -6 erg cm -2 s _I . The spectral slope of vF v cx i/°- 5±0A im- 
plies an X-ray flux of F> s r (20ms) «2x 10 -7 erg cm -2 s _1 in 
the 0.3 - 10 keV XRT range. This, in turn, implies an average 
temporal decay index of (a) w 1.3 — 1.4 between 20 ms and 
200 s. One might expect (a) to be somewhat smaller, as the 
maximal value of a is (3p — 2)/4 (i.e. 1.375 for p = 2.5) 3 at 
v > max(^ m , v c ) which is above 300 keV at 20 ms. This results 
in overproducing the flux at 200 s by a factor of ~ 10-20 for 
p = 2.5. The observed flux is reproduced for p w 2.8. Lower 
values of p might still be possible if, e.g., there are significant 
radiative losses or a much higher £r in the very early after- 
glow. 

4. DISCUSSION AND PROSPECTS 

From the inferred energy per solid angle, simple blast- wave 
models seem able to accommodate the data on the afterglow 
of GRB 050509b. Constraints on the angle-integrated 7-ray 
energy are not very stringent — the outflow could be con- 
centrated in a high Lorentz factor beam only a few degrees 
across, or actually be wider. Standard arguments concerning 
the op acity of a relativistically expanding fireball iPaczvnsQ 
1986) indicate that Lorentz factors T > 10 2 are required, with 
a baryon load ing no larger than ~ IQ^Mq. As we have ar- 
gued in § 13.11 for GRB 050509b internal shocks face the prob- 



3 This is valid before the jet break time 1 Granot & Sari 2002), which most 
likely occurs significantly later than 200 s. 



4 



Lee, Ramirez-Ruiz & Granot 



lem of explaining the observed peak energy. With an external 
shock, the required Lorentz factor is high by usual standards, 
and such accelerations would accordingly require a remark- 
ably low baryon loading close to the central engine. 

Only detailed simulations in full GR will provide us with 
the details of the merger process in a compact binary. How- 
ever, an approximate treatment using variable compressibil- 
ity in the equation of state and a range of mass ratios leads 
to similar outcomes, suggesting that the creation of a dense 
torus is a robust result. If the central engine involves such a 
configuration, is it possible to discriminate between the alter- 
nate modes for its formation: compact merger or collapsar? 
Accurate localizations of further events should help to con- 
firm or reject the latter option, since a collapsar would occur 
in or near a region of recent star formation, contrary to the 
expectations concerning compact object mergers (see § Q. 
A more direct test would obviously be a detect ion, or lack 
there of, of a supe rnova-like signature jBloom et all 12005: 
iHiorth et alJl2005h . Definitive and spectacular confirmation 
could come from the detection of a coincident gravitational 
wave signal in the 0.1-1 kHz range, since mass determina- 
tions in X-ray binaries and the binary pulsars indicate that in 
NS-NS systems mass ratios should be close to unity, whereas 
in BH-NS binaries they should be smaller than 1/3. Accu- 
rate measurement of the inspiral waveform in the LIGO band 
would allow simultaneous determination of the ratio of re- 
duced to total system mass, /i/Mj and of the "chirp" mass 
M c = (miOT2) 3//5 /(Mx) 1 / 5 , from which the mass ratio can be 
derived. 

GRB 050509b is the first event in the short class of bursts 
for which we have an accurate localization and a tentative 
distance indicator, based on an association with an elliptical 

4 It is important to note that the natural time scale for a collapsing envelope 
to produce a GRB is given by the fall-back time, which is longer than a few 



galaxy at z = 0.2248. At the inferred distance of ~ 1 Gpc, we 
have shown here that the energetics and duration can be ac- 
counted for by small, dense disks around stellar mass black 
holes, based on dynamical modeling of such systems. The 
lack of a SN-like signature in the optical at that distance 
(Hiorth et al. 2005) argues against a collapsar/hypernova pro- 
genitor. Putting GRB 050509b at a significantly higher red- 
shift places more serious constraints not only due to the en- 
ergetics, but particularly because of the short duration: at 
Z = 3, -E 7i iso ~ 4.6 x 10 5() erg and tso « 8 ms (see Figures ^ 
and|2j. This is hard to reconcile with the current models, and 
makes it unlikely that a collapsing stellar core is at the ori- 
gin of GRB 050509b. The observed duration distribution of 
bursts may be affected by the mechanism responsible for the 
production of the relativistic outflow, with magnetically pow- 
ered events more faithfully reflecting the intrinsic population. 
GRB 050509b is in many respects an unusual event, being so 
short and apparently sub-energetic. 

Much progress has been made in understanding how 7-rays 
arise from the sudden deposition of energy in a small vol- 
ume, and in deriving the properties of the afterglows that fol- 
low. The identity of short-burst progenitors remains a stand- 
ing mystery, which further observations of events similar to 
GRB 050509b will hopefully help elucidate. 
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